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Abstract 

We formulate and study the nonequilibrium dynamics of strings near the singularity of 
the time-dependent plane wave background in the framework of the Nonequilibrium Thermo 
Field Dynamics (NETFD). In particular, we construct the Hilbert space of the thermal string 
oscillators at nonequilibrium and generalize the NETFD to describe the coordinates of the 
center of mass of the thermal string. The equations of motion of the thermal fields and the 
Hamiltonian are derived. Due to the time-dependence of the oscillator frequencies, a countert- 
erm is present in the Hamiltonian. This counterterm determines the correlation functions in 
a perturbative fashion. We compute the two point correlation function of the thermal string 
at zero order in the power expansion. 
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1 Introduction 



Understanding the dynamics of strings in as large a number of spacetime manifolds as possi- 
ble represents an essential step in the development of the string theory and the design of new 
applications and tests of it. In particular, one expects to obtain important informations con- 
cerning the semiclassical properties of the gravitational interaction in the string theory such 
as the backreaction of the matter to the spacetime background, the cosmological singularity 
and evolution and the initial conditions of the cosmology. Since there are no methods available 
to explore systematically the full string physics in arbitrary spacetimes (see [1] for a review 
of some early results), just a few backgrounds in which the above problems can be addressed 
consistently within the framework of the available semiclassical methods have been studied 
thoroughly. 

The manifolds on which it is possible to approach the string dynamics from the semiclassical 
point of view have the property of possessing time-like Killing vectors at least locally. They can 
be generalized to string backgrounds that are formed by the time-dependent spacetimes and 
massless fields from the closed sector. From these, the distinguished class of time-dependent 
plane wave backgrounds have received a considerably attention lately. There are at least three 
reasons for which the plane waves are interesting for the string theory. In the first place, 
they can be related to cosmological, string and D-brane backgrounds through the Penrose- 
Giiven limit procedure as shown in [2]. Secondly, the plane waves have the interesting property 
that they generate conformal backgrounds in the string theory in which the strings can be 
canonically quantized [SlH]. And thirdly, the techniques developed to quantize the strings in 
plane wave backgrounds can be applied to string theory in other backgrounds with similar 
symmetry structure, for example, to the light-like [5] and the Kasner-like backgrounds [6l [TJ 
|8](see for a recent review [9] and the references therein). The symmetries of the plane wave 
backgrounds allows one to study the problems related to the cosmological singularity in a 
simpler framework in which the tools of the perturbative string quantization are available. 

In order to fully understand the string dynamics near the cosmological singularities or the 
singularities of the time-dependent plane waves, it is imperative to investigate the thermal 
aspects of the string theory in the corresponding time-dependent backgrounds. In the static 
backgrounds, the thermal effects can be studied under the hypothesis of the local thermal 
equilibrium. Consequently, the thermodynamical quantities of interest can be derived from 
the statistical partition function of string. On the other hand, the assumption of the thermal 
equilibrium in the time-dependent plane wave backgrounds leads to the conclusion that the 
string has inequivalent canonical representations near the singularity usually taken at zero or 
infinity time and at an arbitrary finite value of time, respectively |10j . In general, the necessity 
of different theories for different values of the time parameter signals a phase transition in 
between the two values. However, that is not displayed by the thermal equilibrium dynamics. 
Extensive analysis of the temperature under the thermal equilibrium hypothesis in time- 
dependent plane wave backgrounds has been carried out in [TTl [12l [131 El IISl HSl [13 [13 
[19] where it has been shown that what is the equivalent of the Hagedorn temperature in 
the Minkowski space-time could be interpreted either as a limiting temperature or a critical 
temperature or even an oscillating modes dependent temperature. 

The results obtained so far by using different techniques of the quantum fields at finite 
temperature point towards the conclusion that the thermal equilibrium is, at best, a hypothesis 
valid only locally and asymptotically away from the singularity. The very general situation in 
the time-dependent backgrounds is that the interaction between the gravitational field and 
the string involves an exchange of energy at a time dependent rate. As a consequence, the 
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frequencies of the string modes are time-dependent functions of time and the string cannot 
maintain its thermal equihbrium. The time-dependence of oscillator frequencies lead to sensible 
effects in the theory. For example, the correlators of the field theory at finite temperature have 
the following general form 

G{x,y) = Go{x,y) - -signc.(t^ -ty)p{x,y), 

where Go{x,y) is the statistical correlator and C is a closed real-time contour. At thermal 
equilibrium, Go{x,y) ~ p{x,y) while in nonequilibrium Go{x,y) is no longer proportional 
to the spectral function p{x,y). Thus, the heuristic arguments as well as the calculations in 
the equilibrium field theory suggest that, in general, the string ensembles in time-dependent 
backgrounds are nonequilibrium system. 

There are few attempts in the literature to study the nonequilibrium strings even in the 
simpler backgrounds, mainly due to the inner difficulties of the nonequilibrium field theory. 
However, since the strings are solvable by canonical methods in the time-dependent plane wave 
backgrounds one could look for a formulation of strings in terms of nonequilibrium ensembles 
of oscillators. To our knowledge, the few works that address this problem are doing so in 
the framework of one approach, known as the Liouville-von Neumann method |2H [22] (based 
on the generalization of the Lewis- Riesenfeld invariant theorem [23]) has been used recently 
to map the time-dependent quantum mechanics to the Renormalization Group flow and was 
applied to the Penrose-Giiven limit of stacks of D-branes and of the Pilch- Warner solution in 
|24j . The same method was used to analyze the string spectrum in the Penrose-Giiven limit 
of the NS5-brane plane wave background in [25] and to construct the density operator of the 
type IIB superstring in the singular scale-invariant plane wave background in |26j where it was 
shown that the Hagedorn behavior of strings is the same as in the flat spacetime. However, the 
results obtained so far are not complete since they provide a simple description at oscillator 
level without explaining the nonequilibrium dynamics of string fields and the way in which 
the thermal string observables or correlations should be calculated from it. 

In this paper, we propose a new formulation of the nonequilibrium string dynamics in the 
time-dependent plane wave backgrounds based on the Nonequilibrium Thermo Field Dynamics 
(NETFD) [27J. The TFD is a real-time canonical operator formalism in which the statistical 
averages at thermodynamical equilibrium are calculated not by performing the explicit manip- 
ulation of the density operator but rather by averaging the corresponding observables in the 
thermal vacuum that obeys the thermal state condition that is defined on the representation 
space of the TFD (also called thermal space p7|). The mathematical structure of the thermal 
space is that of the direct product of two identical copies of the Hilbert space of the system 
at zero temperature over the field of temperature dependent functions and it is endowed with 
an involution denoted by tilde. The physical interpretation of the thermal vectors is that of 
thermal degrees of freedom expressed in terms of tilde and non-tilde quasi-particles. There is 
a group of linear involution-invariant and time-invariant maps from the direct product of the 
tilde and non-tilde Hilbert spaces at zero temperature to thermal spaces whose elements are 
the thermal Bogoliubov operators that preserve the canonical structure. The thermal group 
is SU{1, 1) for bosons and 5*0(2) for fermions. The choice of a certain Bogoliubov transfor- 
mation is equivalent to the choice of the thermal vacuum state defined by the corresponding 
condensate of pairs of tilde and non-tilde particles. This is so because the Bogoliubov operators 
entangle in a specific way the states from the tilde and non-tilde spaces at zero temperature 
to produces thermal states [27]. The equilibrium TFD was applied to the string theory and 
the string field theory in flat spacetime in[28l[29l[30l|3ll[Ml[331[Ml[35l[Ml[371[^ More 
recently, it has been used to calculate the thermodynamics of strings in static backgrounds in 
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[m 221 251 211 23 23 27] and to construct the thermal D-branes and calculate their thermo- 
dynamical functions in Il8l29l[5Ql[5ll[52l[53l[5l] (see for reviews [551 [Ml EZ] ) . In [58] , the 
equilibrium TFD was generalized to accommodate the topological sector of the strings and 
D-branes on B^'P x T^-p^i. 

The NETFD is a canonical real-time formalism in which the nonequilibrium correlators can 
be calculated from states and operators that have a simple algebraic structure. The NETFD 
formalism retains the basic structure of the equilibrium TFD with the thermal space being 
the direct product of two identical copies of the Hilbert space of the system acted upon by two 
kinds of operators with and without tilde, respectively. They form the representation space 
of the operator algebra. The thermal vectors are interpreted as thermal degrees of freedom of 
the system in nonequilibrium and they can be constructed in the Fock representation of the 
thermal space from the thermal vacuum which is a solution of the master equation 



This is the Schroedinger equation with the infinitesimal time evolution generated by the 
operator [HTl [621163] 



where H and H are the Hamiltonians of the non-tilde and tilde degrees of freedom, respectively. 
The minus sign is due to the unbounded energy of the tilde excitations and is a characteristic 
of the thermal processes in the TFD formalism |27] . Also, it shows that in nonequilibrium 
the thermal vacuum state is generally unstable and time-dependent. Nevertheless, one can 
still define a set of time-dependent creation and annihilation operators such that the thermal 
vacuum be canonically defined for any value of the time parameter |27J. Since the time- 
dependence is carried by the states and the operators, new time-independent operators are 
introduced by time-dependent Bogoliubov maps in order to define stable quasi-particle states. 
They allow one to calculate the correlators of the time-dependent thermal operators in the 
time-independent quasi-particle vacuum with or without interactions [27'. The relativistic 
generalization of the formalism has been proposed very recently in |65t i66j^ . 

This paper is organized as follows. In Section 2 we briefly review the canonical structure 
of the free string in the time-dependent plane wave background following [3]. After gauging 
the world-sheet and the target-space symmetries, the string can be reduced to an ensemble 
of time-dependent non- interacting oscillators in the Brinkmann coordinates. By exploiting 
this structure, we construct in Section 3 the time-dependent Bogoliubov operators and the 
time-independent quasi-particle vacuum. In Section 4 we construct the time-dependent correla- 
tionfunctions of the thermal string. We show that they contain the transition matrix between 
the initial time to = ~oo and the time of the singularity ti = 0, respectively. In NETFD the 
transition matrix is defined in terms of the interaction Hamiltonian of the tilde and non-tilde 
strings and a counterterm that is responsible for the time evolution of the number operators. A 
double Dyson-Schwinger expansion is possible in the two terms. We compute the correlations 
at zero order of the expansion. The last section is devoted to conclusions. In order to make 
the paper self-consistent, some basic facts about the NETFD formalism are collected in the 
Appendix. 

^While we have been finishing this paper, we learned of the reference |67j in which similar results to ours are 
obtained for the complex scalar field in the Minkowski spacetime. The main difference, other than the context, is the 
time-depenence of the frequencies of the string oscillators in contrast to the constant frequencies of the oscillators 
of the scalar field from [67] . 



i^m)) = H\m- 




H = H -H 
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2 Strings in plane wave background 



In this section we review the canonical quantization of the free closed string in the time- 
dependent singular plane wave background [HE] and establish our notations. In what follows, 
we are going to use the Brinkmann coordinates in which the background belongs to the 
following general class 



ds^ = 2dx~^dx 
(p = cl){x~^). 



X{x+) {x'dx+f +{dx'f . 



(3) 
(4) 



Here, x^ are the light-cone coordinates in the target-space and i,j = 2,3, . . . D. The dilaton 
compensates in the Ricci tensor the contribution of the metric and ensures that the background 
is conformally flat. In the case of the singular plane waves, the function A (x) behaves as 
A (x) — ?■ -|- O (x*") at X — ?■ 0, where A; is a real parameter and r > 2. For > the mass 
is positive and the dilaton takes the form 



0(x^ 



(5) 



where C is a constant and D is the spacetime dimension. The string action in curved dilatonic 
background has the general form |68] 



1 



Ana' 



d^ay 



-h 



/i"^ (a) g^, (x) 9„x^ {a) 9^x^ [a] + -R^^^ {a) 4> {a) 



(6) 



where a" = {a^,a^) = (r, a) are the world-sheet coordinates, h'^^ (a) is the world-sheet metric, 
x'^ (a) are the string fields, g^i, (x) is the space-time metric defined by the line element ([3]) with 
= ±,2,3,...!?, and R^"^^ (a) is the scalar curvature of the world-sheet. The symmetries 
of the action ([6]) are fixed by choosing dahaa = and the light-cone gauge x"*" = a'p'^T. 

The general solution of the equations of motion obtained from Q has the Fourier expansion 

m 

x\t, a) 



/ a' \ - 1 r ^ 



n>0 



oi] 2inc! 



(7) 

The creation and annihilation operators from x*(r, cr) are normalized as string modes rather 
than oscillator modes and the following linear superpositions of the Bessel functions have been 
introduced 



ITTV , 



Fnir) = exp( — —)./nTTu 



(27rr) 
(27rr) 



iY„ 



(27rr) 
(27rr) 



(8) 
(9) 



where v = (1/2) (l -|- \/l — Ak) . The modes of the string center of mass can be written in 
terms of time-independent oscillators as follows 



a' 

2{2u-l) 



-l-u 



2i ( a, 



^ \(a 



o + "o 



2i ( — Q, 



Inr 



,k 



(10) 
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The Hamiltonian of string in the hght-cone gauge has the following oscillator structure 



H 



a 

y 



D 

E 

1=2 
D oo 



+ ^ E E ("n < + f^nf^n) " (t) </3; - K (r) ajt^^f 



j=2 n=l 

where the time-dependent coefficients r2„ (r) and <I>n(T) are given by the following relations 



Qn (r) 



1 + 
1 + 



I' 



V 
4^2 



|i^n(T)|V|G„(r)|2 



2nT 



[F„(r)G:(r) + F„*(T)G„(T)] 



2nr 



i^n(r)G„(r). 



(12) 
(13) 



From equations (|12|) and one can see that the string behaves as a collection of self- 

interacting time-dependent string modes even if it is a free string. The non-diagonal inter- 
acting terms do not mix modes of different frequencies. They are generated by the time- 
dependent couplings iln (r) and $n(''") induced by the time-dependent gravitational field. The 
self-interacting sector of the Hamiltonian should be treated non-perturbatively because 0,n (t) 
and ^niT) ai"e of the same strength. However, as was shown in [4J, the Hamiltonian can be 
written as a collection of non-interacting time-dependent harmonic oscillators by applying the 
following linear map 



Ai (r) = n„ (r)< + < (t) (3'^ , 
A^ir) = < iT)atl + Vn (r)/?;, 
BiiT) = Un (t)/3; + < (T)aJ, 
BiHr) = < (r)/3;t + ^;„(r)<, 



where 



Un (r) 
Vn (t) 



2 
1 



2u}n (r) 



drFnir) 



-2iui„{T)T 



-Fnir) + 



2u)n (r) 



The A and B oscillators satisfy the canonical commutation relations 

(r) ,Bil{T)\ = [A'^{T), (^)] = , [yl^ (r ) , 5^ (r)] = [a'^ (r) , i?^t (^) 

In terms of the A and B operators, the light-cone Hamiltonian takes the diagonal form 



(14) 
(15) 
(16) 
(17) 



(18) 
(19) 



= 0. 

(20) 



D 



1=2 



D oo 



+ ^ E E h + ^0] + ^ ' 



i=2 n=l 



where the time-dependent frequencies are 



-|- 



k 
4^' 



(22) 
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The time-dependent function 



hiT) = iD-2)J2^n{r) (23) 



n=l 



plays the role of the normal-ordering constant. It is logarithmic divergent and it can be can- 
celled by a field renormalization of the dilaton for all values of r [4j. The string fields ([7]) can 
be expanded in terms of A and B oscillators as follows 

x\t, a) = xiir) + iMY,^^=r=, {e-'^'^-^^^^Ai (r) e^-^^ - e^^-M^Ajt (r) e-^^-j 
^ ^ V^^n (r) J 

(24) 

+ ij^y^ ^ |e-2^'^"(^)^5; (r) 6-2*"^^ - e2^^"(^)^s;t (r) e^^""^) . 

In general, the particle-like states of time-dependent oscillators are well defined only asymp- 
totically [27]. In the present case, there are two asymptotic Fock spaces of string oscillators 
for observers at r — t- ±00. Since the oscillator sector is asymptotically the same as that of the 
free string in the flat space-time, the asymptotic states can be constructed from the light-cone 
vacua \0{±(X)))^_i^ which are eigenstates of and operators. At finite r, the free string 
picture is lost due to the self- interactions. However, there are free states that can be obtained 
by exciting the time-dependent vacuum |0(r))^_^ in the representation defined by the op- 
erators A and B, respectively. The dynamics of the center-of-mass of string does not reduce 
to the corresponding one in flat space-time, and thus is characteristic to the time-dependent 
plane wave background. 



3 Nonequilibrium dynamics of strings 

The fact that the string theory in the time-dependent plane wave backgrounds is solvable in 
the canonical quantization suggests that the nonequilibrium dynamics could have a canonical 
formulation, too. There are several methods that can be used to this end. As mentioned 
in the introduction, we are going to formulate the nonequilibrium dynamics in the NETFD 
framework because of the simple interpretation that can be given to the states, the operators 
and the nonequilibrium correlators. 

3.1 Thermal string fields 

The first step to be taken in order to develop the NETFD method for strings is to define the 
thermal degrees of freedom. The thermal states belong to the thermal space % = T-L®T-L which 
is the direct product of two identical copies of the Hilbert space of string. On H act operators 
of the form 1 and 1®0 that can be conveniently tensored with vectors from to form 
the so called thermal doublets [2^. We start with the diagonal form of the theory in which the 
string fields are given by the relation (j24p . The coordinates of the tilde string are obtained by 
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applying the tilde operation to x*(r, a) (see the Appendix) 



^ ^ Vn^n (r) J 
;^/Zv^_ _A_ |e2^^"(^)^^; (r) 62*"^^ - e-2*'^"(^)^fl^ (r) 6-2*"^^} . (25) 



n>0 



It is convenient to absorb the time-dependent exponentials into the canonical operators by 
the mean of the following change of variables 

e-^^'^-^^^^Al, (r) = aj, (r) , e^^-M-^^t (^) = ^ (^) ^ (26) 

with similar redefinition of operators in other sectors. The string thermal doublet, denoted by 
(/)*"(r, cr), can be obtained by organizing the fields from (|24|) and (|25|) in to a two dimensional 
vector field in such a way that the a — b structure be maintained 

</.^"(r, a) = C(t) + ^) + ^rir, a) 

= ^h"{r) + ^J^J2^^= {aT (r) e^"^'^ - (s.a'^ (t)^)^-^--} 

- /Z^ |6- (^) - (.36^ (r)^)%2in.| ^ (27) 



n>0 



where a is the index, (^o"(''") stands for the thermal doublet of the center of mass and S3 
is the Pauli matrix. Also, we have introduced the standard thermal doublet 

C(r)= Q^l^^^), C(r) = (a?(r) -<(r)). (28) 

Similar operators can be introduced in the h - sector. The conjugate of the thermal field given 
by the relation ([27j) is obtained by applying the operations defined in the Appendix 



r{r,a) = C(r) + ^(t,^) + 

+ ^^^^Y1 r^-r^ {~^n (r) e^™^ - U {rf ssY e-''^^} . (29) 

Note that the coordinates of the center of mass do not have the standard form of a thermal 
doublet as given by the NETFD. This can be remedied by organizing the coordinates as follows 

4ir) = U{T)al + U*{T)ai\ (30) 

where the time-dependent coefficients are given by the following relations 



U{t)={ V 2(2-1)' 'V2u-1^ , /^^4 _ ^3^^ 

^ - W2a'T\nT, k = l 
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The coordinates of the center of mass of the tilde string can be obtained by applying the 
conjugation rules given in the Appendix. Then the thermal fields are generalized to the center 
of mass by the following relations 



U{T)af + U*{T) {ssaff) 



(32) 



The thermal doublets define a time-parametrized family of representations of the thermal field 
constructed from the thermal vacua {|0(r))}. Any of these representations can be mapped to 
the time-independent quasi-particle representation by the inverse of the Bogoliubov map (see 
the Appendix) 

(33) 



ar (r) 



Then the nonequilibrium dynamics is determined by the time variation of the thermal doublets 
a]^ (r) and b]^ (r) in the time- independent vacuum |0) defined by and Xn^ respectively. The 
Bogoliubov operators can be different in the a and b sectors, respectively, but they are the same 
for all transversal directions. However, since the two sectors differ only in the orientation of 
the modes along the space-like direction of the world-sheet, we can take Ba,n(''") = Bb^„(r) = 
B„(rjl. It follows that the Bogoliubov map of the k - mode has the form given by the relation 
(jl06p in the linear thermal gauge 



a =exp 



= «f (^) exp 



TO 



i I d\LOk{\) 



Xk = exp 



xT = (■^) 



i / d\oJk{X) 



where the positive exponential represents the pure complex phase function and 



nkir)6kd = nl{T)6kd = (0(r) 



0(r) 



(34) 
(35) 
(36) 
(37) 

(38) 



The initial boundary conditions are taken at tq — t- — oo. The time-dependent quasi-particle 
operators are obtained by multiplying the transformations given in the equations (j34[) - (I37[) 
with the inverse of the corresponding phase functions to obtain 



(r) = Bf (r) al^ (r) 
r) 
r) 
r) 



(r), 



Bf (r) 6; 



xT 



r) 

11 /3« 



bfir) [BkT {r 



(39) 
(40) 
(41) 
(42) 



The quasi-particle representation is based on the oscillator equation of motion that is satisfied 
by the string modes ([7]) and by the natural assumption that the system of oscillators evolves 
according to the Schroedinger equation with the Hamiltonian (j2ip in the a — b representation 



^In the case of the complex scalar field discussed in [57] this equality does not hold because the left and right 
moving states have different charges. 
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[Ij. By taking the first derivative of ^^q (r) with respect with r one can easily check that 
a{/3 (r) satisfies the first order differential equation 



dT 



^^kir) +Pk{r) 



al3 



«f (r) = 0, 



(43) 



where 



.dnkij) 
' dr 



(44) 



The same equations hold for bj^f3 (r) operators. This is a general result of NETFD for the set 
of time-dependent oscillators [27] . From (j43p it follows that the time-evolution of the thermal 
string modes obeys the following equations 



, d ja / \ 
(-) 

. d _jQi , , 
(-) 



''k 



a\r [t] ,HQa 



where 



Hgir) = X; E [c^fcCr)^'' + Pf{r)] (r) af (r) . 



(45) 
(46) 

(47) 



i=2 k 



The second term in the equation (j47p is the thermal counterterm which should be added to 
the Hamiltonian as a consequence of the time-dependent Bogoliubov map. As can be seen 
from its definition ()44p it is related to the variation of nfc(r) from the relation (jSSp . The same 
relations holds in the b - sector. Since the total Hamiltonian can be written as 1271 



H = Ho + K 



int 



Hq + Hi, 



(48) 



the counterterm Hq{t) = Hcm{T) + Hq{t) + Hq{t) from (|i8|) should be compensated by 
a term to be added to the interaction term. In the case of the free theory, the interaction 
Hamiltonian contains only the compensator of the counterterm. Note that the quasi-particle 
fields, even if stable, cannot be used to define the asymptotic thermal states because the 
thermal Hamiltonian is unbounded from below [27] . 



3.2 Dynamics of thermal fields 

Once the time-evolution of each mode of the thermal string, described by the equations (j45p 
and (|46p . is understood, we can proceed to deriving the equations of motion of the thermal 
fields. Since we are working in the Hamiltonian formalism, we need to determine the canonical 
conjugates of (/'*"(t, a) and 0*"(t, a). The conjugate momenta are defined such that they satisfy 
the equal-time commutation relations [27j 

(/.^"(r,cj),7r-''^(r, a') 



= i5'^ 5{a - a')5'^P . (50) 



(49) 



<A^-(r,a),7r^-^(r,a') 
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By using the (|45p and (|46|) in the above commutators we obtain 
7r-(r, a) = vr^"(r) + (r, a) + vr^ (r, a) 

= ^[u~\r) (.3a^)" + (C/*(r))-ia- 

+ -L= y/^^^) {a- (r) e^^^ + (.saj, (r)^) ° e'^"-} 

"V^"^ n>0 



-2ina 



n>0 



(51) 



The conjugate momenta of the fields (P^°'{t, a) can be obtained by applying the tilde conjuga- 
tion to TT*" (r, a) and they take the following form 

^-(r, a) = Tti^ir) + a) + 7fr(r, a) 

i 

/2a 



n>0 



^2a' 



n>0 



The fields obtained in this way satisfy the commutation relations 



(52) 



(53) 



which qualifies them for the canonical momenta of the (^*"(t, a) fields. 

Let us determine the equations of motion of the oscillator sector of the fields 
(j)'^°'{T,a),(j)^°'{T,a),TT^°'{T,a) and ■7f*"(r, o"). They can be obtained by deriving the fields with 
respect to r and by using the equations ([^5]) and (06]). Before doing the calculation, it is 
convenient to express the oscillator frequencies in string units through the rescalling 



nUn (r) 



a' 



(54) 



Lengthy calculations give the following equations of motion for the oscillating modes 

"(l + 2i9,n,v|V-iro)"^a, + 5,V- (l + 2ia,n,v|V-iTo)"''] (/.j,i(T, a) = <°,(r, a), 

(1 - 2ia.n,v|V-iro)"^ V2<A^^,(r,a) + [dr - • V-^] vr^ (r,a) = 0. 



(55) 
(56) 



Here, we have introduced the notation V = y —d"^ + The equations ([55]) and ([56|) describe 
the classical dynamics of the oscillator sector of string fields near the plane wave singularity. 
Since the equations of motion are the result of the time evolution generated by the classical 
Hamiltonian, one can use ([5^ . ([^Sj) and to integrate the Hamilton equations. The resulting 
Hamiltonian is composed from two terms 



Hq — Hem + Hose, 



(57) 
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where Hem is the Hamiltonian of the center of mass and 

Hose = ^ daY^Yl (1 + 2z9.n|v|V-iro)°^ <i 

1=2 n>0 

+ <^*°(T,a) (1 - 2i5.n|v|V-iro)°'' VVoi(T,a) 

- (1 + r„)°'' • (1 + 2i5.n|v| V-^To)^^ <A^L(r, cr) 

-<A^"(r,a)a,V-V-Vj,-J, (58) 

where 

A -1\ /2i5,n|v|V-i l-2i9,n|v|V- 

° Vl -ly' " V2i5^n|v|V-i l-2i9^n|v|V" 

The oscillator Hamiltonian can be split in two terms: the free oscillator Hose term and the 
counterterm Q which can be both read off of the right hand side of ()58p . The Hamiltonian 
([571) generates the equations of motion ([551) . (f56]) and (f63]) . The counterterm Q provides the 
dynamics of the number function n|y|(r) from which is derived the time-evolution of each 
oscillator mode. 




3.3 Dynamics of center of mass 

Note that the derivation of the oscillator terms from (j5ip and (j52p follows from the standard 
NETFD method. However, the coefficient of the zero mode component vrg"(r) is not given by 
the general formalism. In order to obtain the correct factor one has to generalize the NETFD 
in the following way. Firstly, we consider the series representation of the delta-function 



1 

2^ 



E 



^2m((T— a) 



(60) 



Then one can show that the commutation relations (|49p and (|50p are reproduced if the coor- 
dinates of the center of mass satisfy the following commutators 



C(T),vr^''(r) 



vr 



From that we can determine the momentum which has the form 



a. 



(61) 



(62) 



This fixes the constant in front of the zero mode term from equation (j5ip . The coefficient of 
the corresponding term from the equation (j52p can be obtained in the same way. 

The equations of motion of the center of mass can be obtained by taking the time derivative 
of the equation ()32p . After some algebra one obtains the following equation 



2tt (drU{T) idrU*{T) 



U*(t) 



"0 



27r - i 27r(27r + i) 
The equation of motion of the momenta is 

1 drU{T) 1 idrU*{T)U{T 



t) = 2ttU{t) 



drU{T) idrU*{T) 

27r — i 2tt + i 



(63) 



Or 



2Tr + i U{t) 



+ 



27r 



[U*{t)Y 



^0 



drU{T) 



1 drU*{T) 



27r(27r + j) [/(r)|C/(r)|2 27r 



i {U*{T)f 
(64) 
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The equations of motion of (^o*(''") ■^o°(''") obtained in the same way. From the 

requirement that the equations of motion be the Hamilton equations for the Hamiltonian of 
the center of mass Hcmt one can calculate the explicit form of it 



D 

E 

i=2 



=.ia-.ta I 



^ ' ^ 2-K-i 2TT + i 

drUir) 1 drU*{T) 



27r (2^ + i) [/(r) I C/(T)r 27r -i{U*{T)y 



1 drU{T) 

27r + i U{t) 
2tt fdrU{T) 



+ 



U*(t) 



1 idrU*{T)U{T) 
2TT-i [C/*(t)]2 

idrU*{T 



j.ia t^\„ia 



27r 



27r(27r + i) 



(65) 



The total Hamiltonian Hem + ^osc determine the dynamics of the thermal string field in 
the proximity of the singularity. 



4 Thermal string correlations 



In order to make predictions about the physical properties of the thermal string, observ- 
able quantities must be defined. At equilibrium, the thermodynamical functions and their 
derivatives provide such observables. In the nonequilibrium, it is more natural to express the 
observables as transition probabilities for string states at different values of the time parame- 
ter. These transition probabilities can be defined in terms of correlations functions. Since we 
are considering the string at zero order in the Euler number, i. e. on cylinder, and in the a — h 
representation there are no interactions among string modes, one would naively expect that 
the correlation numbers be defined by 



(T[0^i"^(ri)0^^-^(T2)---<^*"""(r„)]>. 



(66) 



However, the analysis from the section shows that even without interactions, there is a coun- 
terterm Q that is generated in the Hose as a consequence of the the time evolution of the 
frequencies. In the interacting field theory, Q should be cancelled by self-energy interaction in 
order to leave the Hamiltonian unchanged [27j. Therefore, in order to define the correlation 
functions one has to work in the interaction picture defined by the interaction Hamiltonian 



Hi = Q. 



(67) 



In this representation the oscillators evolve according as given in the equations (I45p and (I46p . 
that is with respect to the free Hamiltonian Hq as defined by the relation ()48p . Then the 
correlation functions are defined as the expectation values 



r[</>*i"i (ri)</.^2"^ (t2) • • • 0*""" (r„)5(0, oo)] 



(68) 



where the S - operator is defined as usual 



S'(0,oo) = lim lim u{ti,T2) = T 

ri— >-0 T2— >oo 



exp {-If d<:,Hi{q] 



(69) 
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The thermal vacuum in the interaction picture satisfies the following relation 



(0|n(0,T) = (0|. (70) 

This is the Schwinger-Dyson formalism for the thermal fields in the canonical formulation. 
The correlations defined by the relation (f68|) can be calculated by expanding in powers of Q 



r[,/.^i"i(n)<^*^"Hr2) • • • (t„)5(0,oo)]^ 

(r[f^i°Hn)</''^"^(r2) • • • 0'"°"(t„)> - i ^r[(/.^i"Hn)</'*^"^(r2) • • • 0'"°"(t„) d^Hi{<,)]^ + ■■■ 

(71) 

Note that the expansion is valid for small values of the time derivatives of Thus, the 
NETFD gives the general method to compute the correlations of the thermal string in a 
perturbative like fashion. 

In order to exemplify this formalism, let us compute the two points correlator at zero 
order. To this end we consider the first term from (j7ip in the time- independent vacuum 



I?*^'"''(Ti,r2;ai,a2) = -i(rU^"(Ti,ai)<^''(T2,a2)lV (72) 



The fields from the equation (l27|l act on the direct space H without mixing the a and h sectors. 
It follows that the Fourier transform of the oscillator propagators has the same form in the 
two sectors and the cross terms vanish. Thus, there are just three non-vanishing terms in the 
left-hand side of the relation ()72p that correspond to the three distinct Hilbert spaces of the 
string fields 

D^'^'^P {n,T2;a,,a2) = D^^^ (n, r2) + Df^^ (n, r2; ai, as) + of^^ (n, rs; ^1,^2) , (73) 

where lower indices denote the center of mass and the corresponding sectors of the string 
oscillators. By substituting the zero modes from (I32p into the above relation, one can show 
that the propagator of the center of mass has the following form 

4''"'' (n, T2) =-\e in - T2) [U{n)U{T2) + U{n)U*{T2)] (I2 + s^r^ 5'^ 

- \0 (t2 - n) [-[/(t2)C/(ti) + U*{T2)U{n)] (I2 - s^r^ (74) 



The calculation of the propagators of the oscillator sectors a and h from the equation (I73p 
is somewhat lengthy but not too difficult. Since the sectors a and h can be related by the 
refiection transformation ui — (T2 — )• CT2 — (Ti one can calculate D^J'°'^ (ri , T2 ; cti , £12 ) and obtain 
jjT-hotP (ri, r2; 0"!, (72) from it by applying the reflection. In order to compute the expectation 
value from -Da '"^ (^i, T2; ci, cs) we use the Fourier decompositions (p7|) and (p9]) and apply 
the Bogoliubov transformations ()34p and ()35p to transform the operators a (r) into the time- 
independent quasi-particle operators ^ that correspond to the vacuum from the equation (j72]) . 
After doing the algebra, the result takes the following form 

n>0 

^g-2^nK-<x,)^J-_ B^{t2)Y^ . (75) 
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We have denoted by the indices + and — terms that would correspond to the advanced and 
retarded propagators, respectively, if the frequency did not depend on time. These terms are 
given by the following equations 



exp 



ird\LOn{\) 



exp -i n dXuJn{X) 



[6 (ri - T2) (I2 + S3)+e (t2 - Ti) (I2 - S3; 



(76) 



2n-s/ujn (n) w„ (T2) 



where we have introduced the matrix for each mode m 

Cmir) -- 



l + 2n^(T) -1 
-1 



(77) 



(78) 



The Fourier transforms of the nonequilibrium propagators in the a sector given by the relations 
(j76p and ()77p are similar to the nonequilibrium propagators of the scalar field in the Minkowski 
spacetime [65l |66l EZ] • As a difference from the scalar field, we note the factor ~\l ujn (ti) a;„ (T2) 
that is not present in the Minkowski spacetime. This factor is a consequence of the dependence 
of the field modes (/'^'^(ti) on ^^JnuJn {t), more specifically,it is the result of the interaction 
between the gravitational background and the strings. The nonequilibrium propagators from 
the b sector can be calculated in exactly the same way. 



5 Conclusions 

In this paper, we have given a canonical method to study the strings near singularity of plane 
wave backgrounds based on the NETFD. This method allows one to construct and interpret 
the Hilbert space and to compute the correlation functions of thermal string. The correlation 
functions of the string oscillators receive corrections from the counterterm Q that generates the 
time evolution of the number parameter n|y|(r). Compared with the literature, these results 
are similar to the ones that are obtained in the case of relativistic scalar field in Minkowski 
[65 11661 [67]. However, the extra degrees of freedom provided by the center of mass of string do 
not have any analogue in the thermal field theory. Therefore, in order to compute their thermal 
behavior, we had to generalize the NETFD by postulating a new commutation relation for 
the canonically conjugate variables of the center of mass ()6ip that is consistent with the equal 
time commutators of the thermal fields given by the relations (j49p . (|50p and (j53p . In the limit 
of small time variations of the number parameter, the NETFD provides a perturbative method 
to compute the correlation function that we have used to determine the two point functions 
at zero order. Due to the importance of the applications of the nonequilibrium string theory, 
the development of the method presented in this paper deserves attention. In particular, it is 
important to calculate the two point functions at first order and the equation that describe 
the evolution of the thermal correlators. We will report on these topics elsewhere. 
Acknowledgements I. V. V. would like to thank to D. Berenstein for discussions and to the 
organizers of LASS 2010 where this work began. Also, I. V. V. would like to thank to J. A. 
Helayel-Neto and A. M. O. de Almeida for hospitality at LAFEX-CBPF where part of this 
work was accomplished. R. N. acknowledges a FAPERJ fellowship. 



14 



A Nonequilibrium Thermo Field Dynamics 

The NETFD is a real-time canonical formalism of thermal quantum field theory. It has 
been originally developed as an alternative to the nonequilibrium imaginary-time formalisms 
which presented difficulties in defining the multi-time correlators and the density operators 
[57]. The NETFD can be derived from the same basic set of axioms as the equilibrium TFD 
which can be summarized as follows 

1. The thermal physical system is described by two sets of commuting (anti-commuting) 
field operators (/>(x) and (pix), respectively, that act on the thermal or total Hilbert 
space which is the direct product of two identical Hilbert spaces 

n = n®n. (79) 

2. To operators O = O (x) , c/)^ (x)) that act on j-L are associated operators O = 
O* (<A(x),0t (^)y 

3. The tilde defines an involution on T-L that obeys the following rules 

(ciOi + C2O2)' = 461 + c^Os, (80) 
(O1O2)' = O1O2, (81) 

'Ot) (82) 



i^Oj =eO, (83) 

for all ci,C2 G C and all 0,0i,02 G End{J-L). e = +1(— 1) for bosonic (fermionic) 
operators. 

4. The vacuum state is invariant under the involution 

\m = \m). m\ = m\- m 

5. The time evolution is generator H should satisfy the following condition 

^iH] =iH. (85) 



From this condition one can see that the Hamiltonian is the difference between the two 
identical Hamiltonians 

H = H -H (86) 
and it is given by the Heisenberg equation for non-tilde operators 

i^p{t)=\0,H\. (87) 
at L J 

6. There is a set of operators f , ^^jthat define the time-independent free quasi- 
particle representation and have the following action on the vacuum 

e|0)=^~|0) = 0, (0|et = (0|et = o. (88) 

The thermal states are defined by the following thermal state condition 

(0| 0{t) = (0| d\t) for bosons, (89) 

(0| 0{t) = (0| d\t) for fermions, (90) 

where Q is determined by the tilde conjugation rules. 
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7. The thermal average of the dynamical observable O is given by 

(O) = (0|O|0). (91) 

8. For any value of the time variable t, there is an invertible map between the total Hilbert 
space and the time-dependent quasi-particle representation given by a time-dependent 
Bogoliubov transformation 

*(«) ^ _ „-u,j ?(«) 



9. The ket-vacuum and the bra-vacuum of the time-dependent representation is defined as 

i{t)m) = Ut)\m) = ^. (93) 

{0(t)|e^(t) = (0(t)||t (t) = o. (94) 
The time-evolution of the time-dependent bra- and ket-vacua is given by the equations 

i||0(t))=^|0(t)), (95) 
(0(t)|i7 = 0. (96) 

10. The stationary thermal states in the Schroedinger representation are defined by the 
equation 

lim ^|0(t)) = 0, (97) 
assuming that the thermal equilibrium is obtained at t — t- oo. 

The NETFD is a canonical formalism. Therefore, there is the time-dependent canonical 
representation defined in terms of the canonical operators {a(t), a^{t), d{t), d^{t)^. The family 
of time-parametrized mappings between the canonical representation and any time-dependent 
quasi-particle representation |^ (t) , (^'^ (t) ,^{t), (i)| is given by the time-dependent invert- 
ible Bogoliubov map 

B{t) : {a(t), at(t), d{t), dt(t)} (t) , (t) , ^ (t) , (t)} . (98) 

The quasi-particle representation is defined by the properties of the Bogoliubov map that 
can satisfy three conditions: i) it can preserve the canonical structure; ii) it can preserve the 
Hermitian conjugation and iii) it can preserve the tilde conjugation: 

B{t){s2(E)s^)B'^{t) = S2®sz, (99) 
B*{t){si®l) = {si(^l)B{t), (100) 
B*{t){siC$si) = {si(^si)B{t), (101) 

where si, 52 and S3 are the Pauli matrices. The first requirement is made in all cases. The 
last two ones define non-standard, albeit useful representations of the thermal field theory. It 
can be shown that the most general form of the non-Hermitian Bogoliubov operator in the 
doublet representation from equation ([92]) is 

B^{t) = (1 + en^ (t))^ e^-W^^ -^i^)'^ ^ (io2) 
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where k denotes the canonical mode of the field, e = 1) for bosons (fermions) and ak, 
nk(t) and 7k(i) are parameters of the gauge group of thermal mappings which is SU{1, 1) for 
bosons and SO{2) for fermions. The pure complex phase function 9k (t) depends on the energy 
of the mode and satisfies the first order differential equation 



. d 



UJn{r)9n (r) . 



The physical interpretation ?ik(i) is that of the number density defined as 



nk(t)5(k-l) = ^0(t) 

The function /^'^ (t) is the statistical distribution 

nk(i) 



0(t) 



1 + enk(t) ' 



(103) 



(104) 



(105) 



The parameters Ok £ [0, 1] lead, in general, to different representations of the thermal field 
in the Liouville space. The choice of the gauge Ok = 1 and 7k (i) = In [1 + enk(t)] renders the 
Bogoliubov transformations linear (|60j) 



BUt) 



1 + e"-k(*) 
—e 



-nk(t) 
1 



(106) 
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